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Deuterated molecules in DM Tau: DCO + , but no HDO 
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Abstract. We report the detection of the J=2-l line of DCO + in the proto-planetary disk of DM Tau and re-analyze the spectrum 
covering the 465 GHz transition of HDO in this source, recently published by Ceccarelli et al. 1 2005). A modelling of the DCO + 
line profile with the source parameters derived from high resolution HCO + observations yields a DCO + /HCO + abundance ratio 
of a 4x 10~ 3 , an order of magnitude smaller than that derived in the low mass cores. The re-analysis of the 465 GHz spectrum, 
using the proper continuum flux (0.5 Jy) and source systemic velocity (6.05 km s _I ), makes it clear that the absorption features 
attributed to HDO and CeH are almost certainly unrelated to these species. We show that the line-to-continuum ratio of an 
absorption line in front of a Keplerian disk can hardly exceed the ratio of the turbulent velocity to the projected rotation velocity 
at the disk edge, unless the line is optically very thick (r > 10 4 ). This ratio is typically 0.1-0.3 in proto-planetary disks and is 
» 0.15 in DM Tau, much smaller than that for the alleged absorption features. We also show that the detection of H2D + in DM 
Tau, previously reported by these authors, is only a 2-sigma detection when the proper velocity is adopted. So far, DCO + is thus 
the only deuterated molecule clearly detected in proto-planetary disks. 
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1. Introduction 

Deuterated molecules are powerful probes of the physical and 
chemical state of star forming regions. They are routinely ob- 
served in low mass pre-stellar cor es where they are f ound to 
be particularly abundant (see e.g. iButner et alJ Il995r) . Since 
the detection of several H-bearing molecules in t he protoplan- 
etary disks surrounding DM Tau and GG Tau JPutrev et alJ 
1997), the question rose of whether deutera ted molecules could 
also b e observed in proto-planetary disks. Ivan Dishoeck et all 
J2003I) have reporte d the detection o f DCO + in the disk of 
TW Hya. Recently, ICeccarelli etaD J2004 120051) (hereafter 
CD04,CD05) have reported the detection in DM Tau of H 2 D + , 
a key ion in deuterium enhancement reactions, and of deuter- 
ated water, HDO. In this letter, we present the detection of 
DCO + in the disk of DM Tau. 

The HDO line reported by CD05 and two lines of CeH 
present in the same spectrum appear in absorption against the 
thermal emission of the dusty disk. These absorptions, if real, 
forces us to revise our views on the disk structure, since wa- 
ter is expected to freeze onto grains at the low temperatures 
(< 25 K) encountered in such disks. They also imply very large 
molecular column densities. We argue in this letter that the low 
signal-to-noise ratio of the observations and an overestimation 
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of the continuum flux by a factor 4 cast doubt on the reality of 
the HDO and CeH detections. 

2. Observational Data 

We have observed the DCO + J=2-l transition at 144.0773 GHz 
using the IRAM 30-m telescope in Oct 199, Feb and Sep 2000, 
using the observing strategy of iDutrev et alJ Jl997l) . The sys- 
tem temperature was about 200 K, and the total integration time 
39 hours. The resulting spectrum is presented in Fig^, with 
the best fit profile derived using a double gaussian with param- 
eters (Vlsr = 5.60 and 6.48 kms" 1 and line width FWHM = 
0.7 kms ' ) taken from the fit of the 13 CO J=2— >1 transition 
given by Guil loteau & Dutrevl Jl994l) . The integrated line flux 
is 0.29 + 0.03 Jy kms" 1 . 

CD05 report the detection of the HDO ln,i — Oo,o line with a 
peak line to continuum ratio of about 0.9, a linewidth of 0.63 ± 
0.18kms _1 , and an LSR velocity of 5.5 kms 1 (using the rest 
frequency of the HDO l ,i -0 ,o transition, 464.9425 GHz from 
|PeLuciaetal]ll97lh . However, the systemic velocity of the 
DM Tau disk is 6.05 + 0.02 km s -1 llGuilloteau & Dutrevl 19 9 3) . 
and all molecular lines detected by Dutr ev et al.l \ 19971) ex- 
hibit a flat-topped or double-peaked profile typical of Keplerian 
disks, with a line width of 1 .5 km s~ . 

Furthermore, CD05 used a continuum flux of 2.0 Jy for DM 
Tau at 464 GHz (650 /mi), based on a model of the source 
presented in their Figure 2. However, the model curve indi- 
cates a 464 GHz flux of * 2 x 10~ 12 ergcrrT 2 s~\ or a 0.5 Jy 
only, in agreement with the two measur ements near t h is fre- 
quency displayed on the same Figure. Darto is etaD d2003l) 
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quote S y (v) = 1 10(v/230GHz) 26 mJy in the millimeter range, 
which yields 630 mJy at 464 GHz. This is an upper limit, since 
the SED is expected to flatten above 250 GHz because of the 
v 2 contribution of the optically thick core (as can be seen in 
Figure 2 of CD051. lAndrews & Williams! ( 2005) measured flux 
densities of 1.08 ± 0.05 Jy at 350 urn and 240 + 10 mJy at 850 
/mi, which, interpolated by a single power law, gives 380 mJy 
(in this case a lower limit to the flux). Thus, the continuum flux 
of DM Tau at 464 GHz is most likely about 0.45 Jy, more than 
a factor of 4 below what has been assumed. FigHJ displays the 
HDO spectrum with the line to continuum ratio computed as- 
suming a continuum flux of 0.5 Jy and an antenna gain of 35 
Jy/K, as appropriate for the JCMT at 460 GHz. 

Using the correct continuum flux and systemic velocity, the 
fitted equivalent width of the HDO line (which is the integrated 
area of the line to continuum ratio) is W — -1.6 + 1.6kms _I . 
Similar results W = -1.2± 1.3 kms -1 are obtained when fixing 
the line width to 1.5 km s 1 , or by using a double -peaked line 
profile as above for DCO + . The absorption feature assigned to 
HDO is thus not significant. 

The inappropriate velocity of 5.5 km s was used by CD04 
for the detection of the ortho-H2D + ly - 1 i o line. When an- 
alyzed with the 6.05 kms~' velocity as above, the integrated 
area of this line is 0.08 : 
tection (see Fig^)). 



0.04 K km s i.e. at best a 2 cr de- 



3. Absorption Lines in Keplerian Disks 

In proto-planetary Keplerian disks, the line formation process 
is different from that in an homogeneous static medium, and 
this for three reasons: 1) the velocity shear created by the disk 
rotation 2) the strong surface density gradient as a function of 
radius, and 3) the radial temperature gradient. 

The first effect has been discussed in deta ils in the con- 
text of line emission from Keplerian disks bv iHorne & M arsh 
(1986), for application to cataclysmic binaries. For absorption, 
the maximum line to continuum ratio is limited by the veloc- 
ity coherence along the line of sight. Unless the disk is seen 
face on, only a fraction of the disk projects at any given ve- 
locity. This fraction can be estimated with simple reasoning. 
From a point at r, 9 (in cylindrical coordinates, with 9 — at 
the intersection of the disk plane with the plane of the sky) in a 
geometrically thin Keplerian disk, the projected velocity along 
the line of sight is 
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Fig. 1. (a) (c) DCO + J=2-l spectrum of DM Tau (this work), 
superimposed on the best fit double-peaked profile, (b) H2D + 
1 i.i - 1 i,o spectrum of DM Tau from CD04, and best fit double- 
peaked profile, (c) HDO spectrum of DM Tau from CD05, 
rescaled with the proper continuum value, superimposed on the 
best fit double-peaked profile (for illustration only, see Sec. 3). 
All fit use the same parameters for velocity and line width. 



Vobs(r, 0) = Vlsr + y/GMJr sin i cos 9 



(1) 




We can set Vlsr = by a simple change of referential. The pro- 
jection of the rotation velocities along the line of sight breaks 
the azimuthal symmetry, but preserves 

V obs (r, n - 6) = -V obs (r, 8) and V obs (r, -9) = V obs (r, 9) (2) 

allowing to study only the interval < 8 < n/2. The locii of 
isovelocity is given by 



r(9) 



(GM,/V 2 bs ) sin 2 ; cos 2 9 



(3) 



Eq|5]implies that the spectrum is symmetric. With a finite local 
line width Av, the line at a given velocity V obs originates from 



Fig. 2. Regions of the disk which yield equal projected veloc- 
ities v b s . The outer ellipse is the projection of the disk outer 
edge. 
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Fig. 3. Predicted absorption line profile as a function of disk parameters. The intensity is represented in units of the continuum 
flux, (a): as a function of outer disk radius (fl ou t 50, 100, 200, 400 and 600 AU). (b) as a function of surface density exponent (p 
— to 3.0 by 0.5) (c): as a function of intrinsic (local, turbulent or thermal) line width (0.05, 0.1, 0.2, 0.4 and 0.8 kms -1 ). (d): 
as a function of molecular column density at 100 AU (10 12 to 10 18 cirT 2 by factors of 10). (e): as a function of disk inclination 
(from thin to thick: 0, 30, 60 and 90°). (f): as a function of excitation temperature (2.7, 3.7, 4.7 and 5.7 K). 



assuming a rectangular line shape for simplification. We define 

vj as : 

v d = VCAf, //? ut sin z (6) 

where fl out is the disk outer radius. In general, unless i is small, 
A v <=; Vd. Then for V ODS < \d + Av/2 , r s (6 = 0) = R out . For 
larger velocities, V b s > \'d + Av/2, r s is smaller and given by 

r s (0 = 0) = GM, sin 2 //(V obs - Av/2) 2 (7) 

Figure|3|indicates the regions of equal projected velocities for 6 
different values: V b s > Vj, Vobs = Vd, Vobs < Vd, and their sym- 
metric counterpart at negative velocities. Fig|2 shows that the 
fraction of the disk covered by the gas at velocities V b s < v d 
is of order Av/v^, and drops very rapidly for larger velocities. 
Accordingly, the line to continuum ratio of absorption lines 
cannot exceed this value, and will be smaller for optically thin 
lines. 

For a more quantitative evaluation, we us ed the radiative 
transfer code developed bv lDutrev et alJ l l 19941) to simulate the 
expected line shape for an absorption line from a Keplerian disk 
with parameters (inclination, dust emission, Av, scale height) 
similar to those of DM Tau. To perform this simulation, we 
assumed that the HDO abundance is constant, and that the 
surface density is a power law as function of radius E(r) = 
SiooC^/lOO A\jy p with exponent p — 1.5. We did not assume 
LTE, but because HDO is difficult to excite by collisions, we 
assumed instead that the levels are populated only by the cos- 
mological background, with a line excitation temperature of 2.7 
K. The continuum emission comes from dust with a power law 
distribution for the temperature T(r) = 20(r/100AU)~ a3 K, 



with the same surface density power law. Unless otherwise 
specified, the disk inclination is set to 30°, the local line width 
to 0.2 km s , the outer radius to 600 AU, the orbital velocity 
at 100 AU to 2.10 kms , corresponding to a stellar mass of 
0.50 Mq, and the HDO column density at 100 AU to 10 14 cm -2 . 

Fig |3^-f show the expected profile for different disk pa- 
rameters. Since the continuum brightness is highly centrally 
peaked, the line to continuum ratio depends significantly nei- 
ther on the disk outer radius, nor on the surface density ex- 
ponent (FigQt-b). Only the separation between the two peaks 
depends on these parameters. Fig|5^ shows that the line shape 
depends strongly on the disk inclination. For a face on disk, the 
line is single-peaked, and its width depends on the local line 
width (with a Vln( T ) opacity broadening factor). For moder- 
ate inclination, the line becomes double-peaked, with the two 
peaks corresponding to the projected rotation velocity near the 
disk outer edge. For an edge-on disk, the opacity towards the 
central star (thus at a zero projected velocity) becomes impor- 
tant, and the central dip is filled again when the column density 
is high enough (as in this example). 

The two most significant effects are due to the intrinsic 
line-width (FigQt) and the molecular column density (Fig|3jl). 
The equivalent width increases nearly linearly with the local 
line width, but only with the logarithm of the column density 
(through the opacity broadening effect), as can be seen in Fig|4] 

We used in Fig|5t-e a minimal excitation temperature, 
2.7 K. For increasing values of T ex , due e.g. to collisions, the 
absorption at zero velocity will first diminish, then reverse to 
emission, the line wing reversing to emission at slightly higher 
values (because the wings originate from the inner disks, where 
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Fig. 4. Left: Equivalent line width of the absorption line as 
function turbulent line width for a column density of 10 x 10 14 
cm 2 . Right: HDO column density for a turbulent line width of 
0.2 km s -1 . Note the logarithmic scale for the column density. 



the dust brightness temperature is higher), as shown in FigFJJ. 
Non uniform mixing with the dust will affect the line intensity, 
althou gh in general weakly, see e.g. IWalterbos & Kennicutll 
dl988h . but will not significantly affect the line shape which 
is dictated by the velocity structure. 

Finally, we note that the simulated profiles could also be ap- 
plied to other molecular transitions, provided the column den- 
sity is changed to the appropriate scale: here, a column density 
of 10 14 cm~ 2 corresponds to a typical opacity of 4 at 100 AU. 

4. Discussion 

We have shown in Sec|5]that there is little evidence that the ab- 
sorption feature identified by CD05 with the loi - Oo,o line of 
HDO (see Fig[T]:) is real and arises in the disk of DM Tau, and 
we have set in Sec|3] further constraints on the expected line 
profiles in such a disk. Two additional absorption features at 
464917.2 MHz and 46505 1 . 1 MHz, observed in the same spec- 
trum, were identified by CD05 with lines of CeH. This identi- 
fication also is almost certainly incorrect. Although two CeH 
transitions are listed at these fr equencies in the JPL Spectral 
line catalog JPickett et alJI 19981) . these are cross-ladder transi- 
tions that connect the 2 Tli/2 upper spin state of C(,H to the 2 Y1t,/2 
ground spin state. These transitions are intrinsically weak and 
originate from levels with high energies (E/k > 80 K): they are 
very unlikely to be seen in DM Tau, especially in absorption. 
Mostly, their frequencies depend on the spin orbit constant A, 
whose value has not been directly measured, a nd are uncertain 
by several hundred MHz (Pears on et al J 19881) . which prevents 
any sensible identification. 

The detection of H 2 D + in DM Tau reported by CD04 is 
also not secure, as it is only at the 2cr level when fitted with 
an appropriate line profile. Since the tentative detection in TW 
Hya reported in the same paper has also not been confirmed, 
but remains consistent with a 2 cr upper limit of 0.3 K kms 1 
from the JCMT (Thi & van Dishoeck 2004, cited by CD04), 
this re-analysis casts doubt on the detection of H 2 D + in proto- 
planetary disks at large. 

This brings the number of deuterated species firmly de- 
tected in proto-planetary disks down to only one, DCO + , in 
two sources: in TW Hya, through the detection of the DCO + 



J=5-4 transition bv lvan Dishoeck et al.ld2003h . and in DM Tau 
(this work). 

We have analyzed the DCO + line using the DM Tau disk 
parameters derived from high angular resolution observations 
of the HCO + J— 1-0 fine performed with the IRAM interfer- 
ometer dPietu et alJl2006l) . We find that this line profile corre- 
sponds to a DCO + /HCO + ratio of 4.0 ± 0.9 x 10~ 3 , provided 
this ratio is constant throughout the disk. Note that a derivation 
which does not take into account the specificity of the line for- 
mation process in a disk (e.g. which would assume HCO + and 
DCO + lines to be optically thin, or would not take into account 
the density and temperature gradients) would give a completely 
different DCO + /HCO + ratio. When no resolved image is avail- 
able, the analysis must still use an approxi mate disk structure, 
as done for example bv lDutrev et"alTdl997t) . The DCO + /HCO + 
abundance ratio derived for DM Tau is smaller than in typical 
dark clouds or (co ld) proto-stellar envelopes (0.01 - 0.07, see 
burner et al.ll995l) . 

Finally, Sec |3] shows that the detection of absorption lines 
from proto-planetary disks require very high sensitivities, and 
that such lines are very insensitive tracers of the total column 
density. At 464 GHz, with 1 kms~' velocity resolution, and a 
system temperature of 400 K, a good sub-mm 12-m antenna 
such as APEX would reach a 1 cr noise of 0.3 Jy in 1 hour. This 
is 3 times more than the maximum absorption signal that could 
be expected in typical proto-planetary disks (line to continuum 
ratio 0.2, flux density 0.5 Jy). A 4 cr detection would require 
more than 100 hours. Only ALMA, with its 50 antennas, would 
reach enough sensitivity, 7 mJy in 1 hour for its minimal 
angular resolution. However, this resolution is already 0.7" at 
this frequency, so that the continuum source would be partially 
resolved. 
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